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We review our recent works on the dynamical localization in the quantum kicked rotator (QKR) 
and the related properties of the classical kicked rotator (the standard map, SM). We introduce the 
Izrailev A-dimensional model of the QKR and analyze the localization properties of the Floquet 
eigenstates [Phys. Rev. E 87, 062905 (2013)], and the statistical properties of the quasienergy 
spectra. We survey normal and anomalous diffusion in the SM, and the related accelerator modes 
[Phys. Rev. E 89, 022905 (2014)]. We analyze the statistical properties [Phys. Rev. E 91,042904 
(2015)] of the localization measure, and show that the reciprocal localization length has an almost 
Gaussian distribution which has a finite variance even in the limit of the infinitely dimensional 
model of the QKR, N —>• oo. This sheds new light on the relation between the QKR and the 
Anderson localization phenomenon in the one-dimensional tight-binding model. It explains the so 
far mysterious strong fluctuations in the scaling properties of the QKR. The reason is that the 
finite bandwidth approximation of the underlying Hamilton dynamical system in the Shepelyansky 
picture [Phys. Rev. Lett. 56, 677 (1986)] does not apply rigorously. These results call for a more 
refined theory of the localization length in the QKR and in similar Floquet systems, where we must 
predict not only the mean value of the inverse of the localization length but also its (Gaussian) 
distribution. We also numerically analyze the related behavior of finite time Lyapunov exponents 
in the SM and of the 2x2 transfer matrix formalism. 
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I. INTRODUCTION 

Quantum chaos, or wave chaos, is the study of the 
phenomena in the quantum domain which correspond to 
the classical chaos in the Hamiltonian systems 01 - Al¬ 
though quantum motion (time evolution of the wavefunc- 
tions) of bound systems with purely discrete energy spec¬ 
trum is ultimately (after a sufficiently long time, asymp¬ 
totically) stable and regular, in fact almost periodic, it 
exhibits many features of the classical motion such as e.g. 
diffusion in a chaotic domain, for time up to the Heisen¬ 
berg time. The Heisenberg time, also called break time, 
is an important time scale in any quantum system, and 
is given by tn = 2 ttH/ AE, where h = 2 ttH is the Planck 
constant and AE is the mean energy level spacing, such 
that the mean energy level density is p(E) = 1/A E. For 
time shorter than approximately tn the quantum diffu¬ 
sion follows the classical chaotic diffusion, but is stopped 
at larger times, just due to the interference phenomena, 
which occur due to the wave nature of the underlying 
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system, and are typically destructive. Pictorially speak¬ 
ing, for time up to tu the quantum system behaves as 
if its evolution operator has a continuous spectrum, like 
the classical one has in the chaotic regime, but at later 
times it senses the discreteness of the spectrum. If the 
quantum diffusion stops, while the classical chaotic dif¬ 
fusion continues, we speak about the dynamical local¬ 
ization, or quantum localization or Chirikov localization , 
first observed in time-dependent systems [3] (see also e.g. 
00)- Through the Fourier transform connection be¬ 
tween the time and energy, the dynamical localization 
reflects itself also in the time-independent eigenfunctions, 
both in the eigenstates of the Floquet operator in time- 
periodic systems and in the eigenfunctions of the time- 
independent, classically chaotic, systems. Namely, if all 
classical transport times like the diffusion time (time nec¬ 
essary to occupy the entire classically available chaotic 
part of the phase space) are all shorter than the Heisen¬ 
berg time tu, we find extended eigenstates, and local¬ 
ized eigenstates otherwise. The subject of this review is 
to summarize our main recent results on dynamical lo¬ 
calization in time-dependent periodic (Floquet) systems, 
exemplified by the quantum kicked rotator (QKR), but 
the approach is nevertheless quite general. 

Another important aspect of quantum chaos is the 
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statistics of the energy spectra of classically chaotic quan- 
tal systems. One of the main cornerstones in the devel¬ 
opment of quantum chaos 0,11 is the finding that in 
classically fully chaotic, ergodic, autonomous Hamilton 
systems with the purely discrete spectrum the fluctua¬ 
tions of the energy spectrum around its mean behavior 
obey the statistical laws described by the Gaussian Ran¬ 
dom Matrix Theory (RMT) H, provided that we are in 
the sufficiently deep semiclassical limit. The latter con¬ 
dition means, as explained above, that all relevant clas¬ 
sical transport times are smaller than Heisenberg time 
tn- This statement is known as the Bohigas -Giannoni 
- Schmit (BGS) conjecture and goes back to their pi¬ 
oneering paper in 1984 7], although some preliminary 
ideas were published in [8j. Since A E oc H? , where / 
is the number of degrees of freedom (= the dimension 
of the configuration space), we see that for sufficiently 
small h the stated condition will always be satisfied. Al¬ 
ternatively, fixing the H, we can go to high energies such 
that the classical transport times become smaller than 
t}j. The role of the antiunitary symmetries that classify 
the statistics in terms of GOE, GUE or GSE (ensembles 
of RMT) has been elucidated in Q, see also 0 , and 
0 , 2 , B, ®. The theoretical foundation for the BGS con¬ 
jecture has been initiated first by Berry llli| . and later 
further developed by Richter and Sieber [12|, arriving fi¬ 
nally in the almost-final proof proposed by the group of 
F. Haake 0®. 

On the other hand, if the system is classically inte- 
grable, Poisson statistics applies, as is well known and 
goes back to the work by Berry and Tabor in 1977 (see 
0, 0 i and the references therein, and for the recent 
advances 0 )- 

In the mixed type regime, where classical regular 
regions coexist in the classical phase space with the 
chaotic regions, being a typical KAM-scenario which 
is the generic situation, the so-called Principle of Uni¬ 
form Semiclassical Condensation (of the Wigner func¬ 
tions of the eigenstates; PUSC) applies, based on the 
ideas by Berry 0, and further extended bv Robnik 
0|. Consequently the Berry-Robnik statistics 0, 0 is 
observed, again under the same semiclassical condition 
stated above requiring that tn is larger than all classical 
transport times. 


The relevant papers dealin g w ith the mixed type 
regime after the work 0 are IBTI-Bsl and the most re¬ 
cent advance was published in [29- 321]. If the couplings 


between the regular eigenstates and chaotic eigenstates 
become important, due to the dynamical tunneling, we 
can use the ensembles of random matrices that capture 
these effects 0]. As the tunneling strengths typically 
decrease exponentially with the inverse effective Planck 
constant, they rapidly disappear with increasing energy, 
or by decreasing the value of the Planck constant. In such 
case the regular and chaotic eigenstates can be separated 
and the dynamical localization in the chaotic eigenstates 
can be studied. For an excellent review of dynamical 
localization in the time-independent systems see the pa¬ 


per by Prosen [0] and the references therein. In such 
a situation it turns out that the Wigner functions of the 
chaotic eigenstates no longer uniformly occupy the entire 
classically accessible chaotic region in the classical phase 
space, but are localized on a proper subset of it. Indeed, 
his has been analyzed with unprecedented precision and 
statistical significance by Batistic and Robnik [Bill. liiTl. [32) 
in case of mixed type systems. The important discovery 
is that the level spacing distribution of the dynamically 
localized chaotic eigenstates in periodic as well as time- 
independent systems is exceedingly well described by the 
Brody distribution, introduced in J35;| , see also 0, with 
the Brody parameter values /3br within the interval [0,1], 
where /3 br = 0 yields the Poisson distribution in case of 
the strongest localization, and /3br = 1 gives the Wigner 
surmise (2D GOE, as an excellent approximation of the 
infinite dimensional GOE). The Brody distribution was 
found to fit the empirical data much better than the dis¬ 
tribution function proposed by Izrailev (see [37, 38] and 
the references therein) characterized by the parameter 
(3iz- It is well known that Brody distribution so far has 
no theoretical foundation, but our empirical results show 
that we have to consider it seriously thereby being moti¬ 
vated for seeking its physical foundation. 

In this review of our papers [398 - 140 we explore the 
quantum kicked rotator (QKR) introduced by Casati et 
al. ® from the classical point of view (the standard 
map, SM), analyze the quantum analog using the N- 
dimensional model of Izrailev, and consider the semiclas¬ 
sical connection between the two pictures. We shall treat 
the cases with the classical dimensionless kick parameter 
K in the range K £ [5,35], and for some purposes even 
up to K = 70, and in the end shall focus on the case 
K = 10, which is the most chaotic one in the sense that 
it is fully chaotic with minimal (in fact undetected) reg¬ 
ular regions among all cases I\ £ [5, 70] and among them 
best exhibits normal diffusion. Izrailev’s IV-dimensional 
model introduced and discussed in 0, [43-140 is treated 
for various N < 3000, which in the limit N —> oo tends 
to the QKR. Due to the finiteness of N the observed 
(dimensionless) localization length of the eigenfunctions 
in the space of the angular momentum quantum number 
does not possess a sharply defined value, but has a cer¬ 
tain distribution instead. Its reciprocal value is almost 
Gaussian distributed. This might be expected on the 
analogy with the finite time Lyapunov exponents in the 
Hamiltonian dynamical systems. In order to corroborate 
the theoretical findings on this topics we perform the nu¬ 
merical analysis of the finite time Lyapunov exponents 
in the standard map (classical kicked rotator), especially 
the decay_of the variance. Indeed, in the Shepelyansky 
picture [0 the localization length can be obtained as the 
inverse of the smallest positive Lyapunov exponent of a 
finite 2fc-dimensional Hamilton system associated with 
the band matrix representation of the QKR, where k is 
the quantum kick parameter (to be precisely defined be¬ 
low). In this picture, N plays the role of time. However, 
unlike the chaotic classical maps or products of transfer 
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matrices in the Anderson tight-binding approximation, 
where the mean value of the finite time Lyapunov ex¬ 
ponents is usually equal to their asymptotical value of 
infinite time and the variance decreases inversely with 
time, as we also carefully checked, here the distribution 
is found to be independent of N: It has a nonzero vari¬ 
ance even in the limit N —> oo. The reason is that the 
quantum kicked rotator at N = oo cannot be exactly 
modeled with finite bandwidth (equal to 2k) band ma¬ 
trices, but only approximately, such that the underlying 
Hamilton system of the Shepelyansky picture has a grow¬ 
ing dimension with N, implying asymptotically an infi¬ 
nite set of Lyapunov exponents and behavior different 
from the finite dimensional Hamiltonian systems. The 
observation of the distribution of the localization length 
around its mean value with finite variance also explains 
the strong fluctuations in the scaling laws of the kicked 
rotator, such as e.g. the entropy localization measure as 
a function of the theoretical scaling parameter A, to be 
discussed below. On the other hand, the two different 
empirical localization measures, namely the mean local¬ 
ization length as extracted directly from the exponen¬ 
tially localized eigenfunctions and the measure based on 
the information entropy of the eigenstates, are perfectly 
well linearly connected and thus equivalent. Therefore 
these results call for a refined theory of the localization 
length in the quantum kicked rotator and similar sys¬ 
tems, where we must predict not only the mean value 
of the inverse localization length but also its (Gaussian) 
distribution, in particular the variance IhI (Manos and 
Robnik 2015). 

The paper is organized as follows. In Sec. [TT] we in¬ 
troduce and study the the classical kicked rotator (stan¬ 
dard map, SM) and concentrate on the role of acceler¬ 
ator modes for the anomalous diffusion. In Sec. uni we 
introduce the quantum kicked rotator (QKR) and the 
A-dimensional Izrailev model and study the localization 
properties, including the scaling laws. In Sec. II VI we an¬ 
alyze the numerical results showing that the localization 
measure has a distribution whose variance does not go to 
zero in the limit N —> oo. In Sec. 0 we discuss the re¬ 
sults and draw the main conclusions, pointing out the 
important differences between the QKR and the one¬ 
dimensional Anderson localization in the tight-binding 
approximation. 

II. THE CLASSICAL KICKED ROTATOR: THE 
STANDARD MAP 

The kicked rotator was introduced by Casati et al. Q 
and is one of the key model systems in classical and quan¬ 
tum chaos, especially for time-periodic (Floquet) sys¬ 
tems. The Hamiltonian function is 

H = + V 0 S T (t) cos 9. (1) 

Here p is the (angular) momentum, I the moment of in¬ 
ertia, Vq is the strength of the periodic kicking, 9 is the 


(canonically conjugate, rotation) angle, and Srit) is the 
periodic Dirac delta function with period T. Between the 
kicks the rotation is free, therefore the Hamilton equa¬ 
tions of motion can be immediately integrated, and thus 
the dynamics can be reduced to the standard map (SM), 
or so-called Chirikov-Taylor map, given by 

Pn+i =Pn + V 0 sin 9 n+ 1 , 

T (2) 

9n +1 = 9 n A —Pm 

and introduced in [Hj-|48j]. Here the quantities (9 n ,p n ) 
refer to their values just immediately after the n-th 
kick. Then, by introducing new dimensionless momen¬ 
tum P n = p n T / 1, we get 

j P n +i = Pn A K sin 0 n .y. i, 

\ $ra+l = On A P n , 

where the system is now governed by a single classical 
dimensionless kick parameter I\ = VoT/I , and the map¬ 
ping is area preserving. 

The generalized diffusion process of the standard map 
[Eq. ©] is defined by 

<(AP) a > =£>„(*>", (4) 

where n is the number of iterations (kicks), and the ex¬ 
ponent p, is in the interval [0,2), and all variables P, 9 
and K are dimensionless. Here DffK) is the generalized 
classical diffusion constant. In the case p = 1 we have 
the normal diffusion , and D\{K) is then the normal dif¬ 
fusion constant, whilst in the case of anomalous diffusion 
we observe subdiffusion when 0 < p < 1 or superdiffu¬ 
sion if 1 < p < 2. In the case p = 2 we have the ballistic 
transport which is associated strictly with the presence 
of accelerator modes. 

In the case of the normal diffusion p = 1 the theoretical 
value of Di(K) is given in the literature, e.g. in [38| or 
0 , 

D (K] = { \ K2 t 1 - 2 UK) (1 - UK))\ , if K > 4.5 
U 10.15(AT - K cr ) 3 , if K cr < K < 4.5 

(5) 

where K cr ~ 0.9716 and J 2 i.IV) is the Bessel function. 
Here we neglect higher terms of order K~ 2 . However, 
there are many important subtle details in the classical 
diffusion further discussed below. 

The dependence of the diffusion constant for the 
growth of the variance of the momentum on AT is very 
sensitive, and described in the theoretical result [Eq. j5j)], 
and fails around the period 1 accelerator mode intervals 

(27m) < K < yj (27rn) 2 A 16, (6) 

n any positive integer. In these intervals for the acceler¬ 
ator modes n = 1 we have two stable fixed points located 
at p = 0, 9 = 7T — 9o and p = 0, 9 = n + 9o, where 
9q = arcsin(27r/A'). There are two unstable fixed points 
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FIG. 1. (Color online) The classical diffusion constant D e g = 
((A P) 2 )/n for the standard map as a function of K (SK = 
0.05) for three discrete times n, i.e., the number of the itera¬ 
tions of the standard map, n = 1000 (lower red dashed line), 
n = 5000 (intermediate blue solid) and n = 10000 (upper 
black dot-dashed). In the background we have plotted the 
classical diffusion constant Di (gray dotted line) [Eq. ([5])]. 
The presence of accelerator modes at certain intervals of K 
(and the sticky objects around) generate anomalous diffusion 
which is rendered by peaks. Here we used sslOOOOO (314x314) 
initial conditions uniformly distributed in a grid on the entire 
phase space [0, 2n] x [0, 2^ r]. 


at p = 0, 9 = 9 o and p = 0, 9 = 2ir — 9q. For example, 
in the case K = 6.5 we have 9q ss 1.31179. Moreover, as 
the diffusion might even be anomalous, we have recalcu¬ 
lated the effective diffusion constant D e g = ((A P) 2 )/n 
numerically, which in general is not equal to the D^ de¬ 
fined in Eq. ©. In Fig. □ we show the D e g for the stan¬ 
dard map as a function of I\ for three discrete times 
?z, i.e., the number of the iterations of the standard map, 
n = 1000 (lower red dashed line), n = 5000 (intermediate 
blue solid) and n = 10000 (upper black dot-dashed). In 
the background we have plotted the theoretical diffusion 
constant Di taking into account only the normal diffusion 
(gray dotted line) [Eq. (0]. The presence of accelerator 
modes at certain intervals of K (and the sticky objects 
around) generates anomalous diffusion which is rendered 
by peaks. Here we used ~ 100000 (314 x 314) initial con¬ 
ditions uniformly distributed in a grid on the entire phase 
space [0,27 t] x [0, 2tt). We see that the dotted theoreti¬ 
cal curve stemming from Eq. © describes the diffusion 
constant well outside the accelerator mode intervals. In 
general, however, the diffusion might be non-normal, de¬ 
scribed in Eq. ©. There are also accelerator modes of 
higher period (2,3,4...) observed and examined below. 

In Fig. [5] we show the variance of the momentum P in 
the standard map [Eq. ©] with K = 6.5 (red crosses) 
where small islands and accelerator mode of period 1 are 
present and I\ = 10.0 (blue stars) where the phase space 
is fully chaotic for the same initial conditions as in Fig. [1] 
as a function of the discrete time n (number of iterations), 



FIG. 2. (Color online) The variance of the momentum P in 
the standard map [Eq. ©] with K = 6.5 (red crosses) where 
small islands and accelerator mode of period 1 are present and 
K = 10.0 (blue stars) where the phase space is fully chaotic for 
the same initial conditions as in Fig.[Qas a function of the dis¬ 
crete time n (number of iterations), in log-log representation. 
The two slopes associated with different types of diffusion are 
p(K = 6.5) = 1.61252 (dotted), p{K = 10.0) = 0.991334 
(solid) with standard deviation errors ±0.01271 (0.7881%) 
and ±0.0009537 (0.0962%) respectively. 


in log-log representation. The two slopes associated with 
different types of diffusion are p(K = 6.5) = 1.61252 
(dotted), p(I\ = 10.0) = 0.991334 (solid) with standard 
deviation errors ±0.01271 (0.7881%) and ±0.0009537 
(0.0962%) respectively. In this manner we have calcu¬ 
lated the diffusion exponent p for all K on the interval 
K e [ K cr , 70] and the result is shown in Fig. [3] We show 
the diffusion exponent p as a function of K after n = 5000 
iterations, using a fine grid of 314 x 314 (« 100000) ini¬ 
tial conditions on the plane (9,P) = (0,2n). The p is 
calculated by the slopes, of the lines of the variance of 
the momentum P as a function of iterations, as it is de¬ 
scribed in [Eq. ©] and for a grid of cells on the entire 
phase space. The intervals on the black horizontal line 
p = 0.9 indicate the intervals of stable accelerator modes 
of period 1 [Eq. ©]. All intervals of K with exponent 
p ss 1 are associated with normal diffusion processes. 
The large peaks (appearing mainly for K > 27t marked 
with full black circles) reflect the anomalous diffusion due 
to accelerator modes [of period 1, being located inside the 
intervals predicted by the Eq. ©]. However, there is a 
number of relatively smaller peaks for K < 27t (more 
clearly presented in the inset panel of Fig. El), whose ori¬ 
gin is accelerator modes of higher period as we will see 
later, and also for 27 t < K < 47r, both these sets are 
marked with empty circle. With the symbol (x) we mark 
few typical examples, close to those peaks, for which the 
diffusion is normal and are also studied in detail in this 
section. 

All the large peaks for K > 2ir , marked with full black 
circles in Fig. [3l correspond to regimes with accelerator 
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modes of period 1 and they decrease monotonically as a 
power law 

f{x) = ax b , (7) 

where a = 2.41645 and b = —0.195896 [blue dotted 
line in Fig. [3] with asymptotic standard error ±0.04294 
(1.777%) and ±0.00537 (2.741%) respectively] indicat¬ 
ing that for K > 70 their effect decreases significantly. 
On the other hand, the size of the successive accelera¬ 
tor modes of period 1 intervals decays with a power law 
defined simply and analytically by the Eq. (j6|) . In or¬ 
der to understand the effect of the presence of acceler¬ 
ator modes in the diffusion and transport properties of 
the phase space in the standard map, we first picked an, 
as much as possible, representative sample of A'-values. 
In more detail, we included in our test-cases all the K- 
values which correspond to all the main peaks appearing 
in Fig. [3] with I\ > 27t together with a few cases from 
the ‘plateaus’ of this curve. Furthermore, we took into 
account the peaks occurring for 1 < K < 2n (see the 
empty black circles in the inset zoom in Fig. E| which 
are associated with accelerator modes of higher periodic¬ 
ity, as it will be seen thereupon. The case with K = 3.8, 
whose /i value is ss 1, is chosen for comparison reasons 
from the plateau and as it turns out has no accelerator 
modes in its phase space causing anomalous diffusion. 
Here we should stress that we repeated the same proce¬ 
dure for larger number of iterations n = 10000 and it 
turns out that the exponent /i has well converged to the 
values shown in Fig. El 

The distribution of the momenta in the case of nor¬ 
mal diffusion is found to be perfect Gaussian, whilst for 
anomalous diffusion a strong departure from the Gaus¬ 
sian distribution is observed, being well fitted by a sta¬ 
ble Levy distribution, characterized by the parameter 
a £ [0,2]. The details are given in reference [40j |. For 
each one of the K -values of the nonlinearity kick param¬ 
eter of Fig. El we have performed a thorough study by 
calculating and comparing the following quantities 

(a) The index of stability a-parameter of the Levy sta¬ 
ble distribution. 

(b) The diffusion exponent p as described in Eq. f4|. 

In the case of normal diffusion (Gaussian statistics) 

for the above quantities, one expects to find a = 2 
for the Levy stable distribution and diffusion exponent 
/i = 1, while in the general case we find other values. We 
have calculated the /i exponent emerging from a small 
box/ensemble of initial conditions in the phase space 
and thus produced the /r-landscape in the phase space of 
the SM. Furthermore, we have also employed the GALI- 
method [E3, IHU to calculate the GALI-index in order to 
identify the regular and chaotic regions in the phase space 
of the SM, which also quantifies the degree of chaos (in¬ 
directly). The GALI-landscape and the /x-landscape are 
found to correspond very well to each other, where it 
seems that the ^-plot contains more information than 
GALI-plot. 


The conclusion of this section is that the SM exhibits 
normal diffusion for most of the K values on the inter¬ 
val [K cr , 70], except for the accelerator mode intervals 
where anomalous diffusion is observed, with the expo¬ 
nent n typically being larger than 1. Using these plots 
and the described methodology we found that the case 
K = 10 is the closest to full chaos (no regular islands 
present) and exhibits the normal diffusion for all initial 
conditions in the phase space of the SM. 


III. THE QUANTUM KICKED ROTATOR AND 
THE IZRAILEV MODEL 


The quantum kicked rotator (QKR) is the quantized 
version of Eq. 0. namely 


h 2 d 2 

H = ~2id(p +v ° 6t ^ cose - 


( 8 ) 


The Floquet operator F acting on the wavefunctions 
(probability amplitudes) 6 € [0, 2n), upon each pe¬ 

riod (of length T) can be written as (see e.g. Iljj, Chapter 

4) 


a ( iVo fl \ ( ihT d 2 \ 

F = exp cos V exp \~~2TW) ’ 


(9) 


where now we have two dimensionless quantum control 
parameters 


k = 


Vo 
h 1 



( 10 ) 


which satisfy the relationship K = kr = VqT/I , K being 
the classical dimensionless control parameter of Eq. ©• 
By using the angular momentum eigenfunctions 


I n) = a n {9) 


-== exjp(inO), 


(ii) 


where n is any integer, we find the matrix elements of F, 
namely 


F„ 


(m|F|n) = exp 



1 T (k\ 

°n—m v'V i 


( 12 ) 


where J l/ (k) is the r'-th order Bessel function. For a 
wavefunction ip(0) we shall denote its angular momen¬ 
tum component (Fourier component) by 

U n = (n\ip) = / a* (d)ip(6) do = 
o 

ift(9) exp (—inO) dO. (13) 

The QKR has very complex dynamics and spectral prop¬ 
erties. As the phase space is infinite (cylinder), p £ 
(—oo, ±oo ),0 £ [0, 27t), the spectrum of the eigenphases 
of F, denoted by (j) n , or the associated quasienergies 
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FIG. 3. (Color online) The diffusion exponent p as a function of K after n = 5000 iterations and for «100000 (314 x 314) 
initial conditions on the plane (9, P) = (0, 2ir). The intervals on the black horizontal line p = 0.9 indicate the intervals of stable 
accelerator modes of period 1 [Eq. ©]• All intervals of K with exponent p sa 1 are associated with normal diffusion processes. 
The large peaks (appearing mainly for K > 2n marked with full black circles) reflect the anomalous diffusion accelerator modes 
(mainly of period 1). The smaller peaks for K < 2n (more clearly presented in the inset panel) originate by accelerator modes 
of higher period together with those for 2n < K < marked with empty circle and a few typical examples close to those 
peaks [marked with the symbol (x)], for which the diffusion is normal, are studied thoroughly later on. The blue dotted line 
corresponds to the power law which describes the decay of the exponent p of the main peaks’ amplitude due to accelerator 
modes of period 1 (see text for more details). 


fujj n = h<j> n /T, introduced by Zeldovich ]52|, can be con¬ 
tinuous, or discrete. It is quite well understood that for 
the resonant values of r 


Airr 



(14) 


q and r being positive integers without common fac¬ 
tor, the spectrum is continuous, as rigorously proven by 


Izrailev and Sliepelyansky [53H56( |. and the dynamics is 
(asymptotically) ballistic, meaning that starting from an 
arbitrary initial state the mean value of the momentum 
( p ) increases linearly in time, and the energy of the sys¬ 
tem E = ( p 2 )/(2I ) grows quadratically without limits. 
For the special case q = r = 1 this can be shown elemen¬ 
tary. Such behavior is a purely quantum effect, called 
the quantum resonance. Also, the regime of quadratic 
energy growth manifests itself only after very large time, 
which grows very fast with the value of the integer q from 
Eq. ([HI) , such that for larger q this regime practically 
cannot be observed. 

For generic values of t/(47t), being irrational number, 
the spectrum is expected to be discrete but infinite. But 
the picture is very complicated. Casati and Guarneri 


57] have proven that for t/(4tt) sufficiently close to a 


rational number, there exists a continuous component 
in the quasienergy spectrum. So, the absence of dy¬ 
namical localization for such cases is expected as well. 


Without a rigorous proof, we finally believe that for all 
other (“good”) irrational values of r/(47r) we indeed have 
discrete spectrum and quantum dynamical localization. 
In such case the quantum dynamics is almost periodic, 
and because of the effective finiteness of the relevant set 
of components u n and of the basis functions involved, 
just due to the exponential localization (see below), it is 
even effectively quasiperiodic (effectively there is a finite 
number of frequencies), and any initial state returns after 
some recurrence time arbitrarily close to the initial state. 
Thus the energy cannot grow indefinitely. 

The asymptotic localized eigenstates are exponentially 
localized. The (dimensionless) theoretical localization 
length in the space of the angular momentum quantum 
numbers is given below, and is equal (after introducing 
some numerical correction factor a M ) to the dimension¬ 
less localization time fw [Eq. (ITT)]) ]. We denote it unlike 
in reference fH] and [39] by C. Therefore, an exponen¬ 
tially localized eigenfunction centered at m in the angular 
momentum space [Eq. dill) ] has the following form 


z exp 



(15) 


where u n is the probability amplitude [Eq. (1131) ] of the 
localized wavefunction ip(0). The argument leading to 
fi oc in Eq. (1161) originates from the observation of the dy¬ 
namical localization by Casati et al. [3j, and in particular 
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from [58j |. and is well explained in [li], in case of normal 
diffusion /z = 1, whilst for general p we gave a theoreti¬ 
cal argument in [39]. We shall denote a = 2/£, and will 
later on determine the a’s directly from the individual 
numerically calculated eigenstate. 

The question arises, where do we see the phenom¬ 
ena (spectral statistics, namely Brody-like level spacing 
distribution) analogous in the quantum chaos of time- 
independent bound systems with discrete spectrum? To 
see these effects the system must have effectively finite 
dimension, because in the infinite dimensional case we 
simply observe Poissonian statistics. Truncation of the 
infinite matrix F mn in Eq. m in tour de force is not 
acceptable, even in the technical case of numerical com¬ 
putations, since after truncation the Floquet operator is 
no longer unitary. 

The only way to obtain a quantum system which shall 
in this sense correspond to the classical dynamical sys¬ 
tem [Eqs. (ID) . OT) and ©] is to introduce a finite N- 
dimensional matrix, which is symmetric unitary, and 
which in the limit N — > oo becomes the infinite dimen¬ 
sional system with the Floquet operator [Eq. (|9jl]. The 
semiclassical limit is k —> oo and r —> 0, such that 
K = kr = constant. As it is well known j38|, for the 
reasons discussed above, the system behaves very simi¬ 
larly for rational and irrational values of t/(47t). Such a 
TV-dimensional model [37] will be introduced below. 

Following [3j| [58} we find that the dimensionless 
Heisenberg time, also called break time or localization 
time, denoted by fi oc , in units of kicking period T, is 
equal to the dimensionless localization length C as shown 
in Eq. (l20l) . 


C 



(16) 


where a M is a numerical constant to be determined em¬ 
pirically, and in case of normal diffusion /z = 1 is close to 
1/2. Since this semiclassical approach and derivation is 
quite important, we repeat the arguments given in [3f|. 

The generalized diffusion process of the standard map 
[Eq. ©] is defined by Eq. ©. As the real physical angu¬ 
lar momentum p and P are connected by P = pT //we 
have for the variance of p the following equation 

((A P) 2 ) = ^D (17) 

(Please note that in ref. [39| there is a typing error in 
formula (13), where T and I must be switched.) Now 
we argue as follows: The general wisdom (golden rule) in 
quantum chaos is that the quantum diffusion follows the 
classical diffusion up to the Heisenberg time (or break 
time, or localization time), defined as 


tH = 


27 tH 
~AE’ 


(18) 


where A E is the mean energy level spacing. In our 
case we have the quasienergies and A E = HAuj , where 


Aw = A<//T, and A (f> is the mean spacing of the eigen- 
phases. This might be estimated at the first sight as 
A(p = 27 t/N, but this is an underestimate, as effectively 
we shall have due to the localization only C levels on the 
interval [0, 27r). Therefore Acf = 27 t/C and we find 

tH = ^ =TC - (19) 

Since T is the period of kicking, and tn is the real phys¬ 
ical continuous time, we get the result that the discrete 
time [number of iterations of Eq. © at which the quan¬ 
tum diffusion stops], the localization time t\ oc is indeed 
equal to the localization length in momentum space, i.e. 

tioc « C. (20) 

Since our derivation is not rigorous, we use the approx¬ 
imation symbol rather than equality, in particular as 
the definition depends linearly on the definition of the 
Heisenberg time. Now the final step: By inspection of the 
dynamics of the Floquet quantal system [Eqs. ©,©] one 
can see (see also the derivation in the Stockmann’s book 
Q) that the value of the variance of the angular mo¬ 
mentum at the point of stopping the diffusion ((A p) 2 ) 
is proportional to H 2 C 2 , and to achieve equality we in¬ 
troduce a dimensionless numerical (empirical,correction) 
factor by writing ((A p) 2 ) = h 2 C 2 which on the 
other hand must be equal just to the classical value at 
stopping time t\ oc , namely equal to (I/T) 2 D^C^. From 
this it follows immediately Eq. m- The numerical con¬ 
stant a M is found empirically by numerical calculations, 
for instance in the literature the case K = 5 with p = 1 
is found to be oq = 0.5 (however, we find numerically 
oq = 0.45, taking into account Eq. (TT6l) when studying 
the model’s localization properties). Thus, we have the 
theoretical formula for the localization length in the case 
of generalized classical diffusion [Eqs. (|H). (fT7|) ]. which we 
use in defining the scaling parameter A below. 

The motion of the QKR [Eq. ©] after one period T of 
the if wavefunction can be described also by the following 
symmetrized Floquet mapping, describing the evolution 
of the kicked rotator from the middle of a free rotation 
over a kick to the middle of the next free rotation, as 
follows 

if{6,t+T) = Uif{6,t), (21) 

f T (Thd 2 \ ( .Ho A (Thd 2 \ 

u = exv \ l uW 2 ) exp C0S V exp VJTW) ■ 

Thus, the ip(0,t) function is determined in the middle of 
the rotation, between two successive kicks. The evolu¬ 
tion operator U of the system corresponds to one period. 
Due to the instant action of the kick, this evolution can be 
written as the product of three non-commuting unitary 
operators, the first and third of which correspond to the 
free rotation during half a period G(r/2) = exp , 

r = hT/I , while the second B{k) = exp (—ik cos 9), 




k = Vq/H, describes the kick. Like before, we have 
only two dimensionless parameters, namely r and k, and 
K = kr = VqT/I. In the case K = fcr 1 the motion is 
well known to be strongly chaotic, for K = 10 certainly 
without any regular islands of stability, as mentioned, 
and also there are no accelerator modes, so that the dif¬ 
fusion is normal (p = 1). We have carefully checked that 
the case K = 10 is the closest to the normal diffusion 
/i = 1 for all K £ [K C r, 70]. The transition to classi¬ 
cal mechanics is described by the limit k —> oo, r —> 0 
while K = const. We shall consider the regimes on the 
interval 3 < k < 20, but will concentrate mostly on the 
semiclassical regime k > K, where r < 1. 

In order to study how the localization affects the sta¬ 
tistical properties of the quasienergy spectra, we use the 
model’s representation in the momentum space with a 
finite number TV of levels [§3: H3, l43-[i3| , which we refer 
to as Izrailev model 

N 

Unit + T) = ^ UnmU m (t), U, TO = 1, 2, ..., TV . (22) 

771=1 


The finite symmetric unitary matrix U nm determines the 
evolution of an TV-dimensional vector, namely the Fourier 
transform u n {t) of ip(9, t), and is composed in the follow¬ 
ing way 


Un.m. 


n'm' 


Gri.m.' Bri/rn/ G 


n'm' 'GT'n'm 5 


(23) 


where Gw = exp (irl 2 / 4) 5w is a diagonal matrix corre¬ 
sponding to free rotation during a half period T / 2, and 
the matrix B n i m i describing the one kick has the follow¬ 
ing form 


Bn 


1 


2TV+ 1 


2JV+1 

E 

i=i 


x exp 


cos 


(n' — to') 


27rZ 


—ik cos 


2TV + 1 

27 rl 


— cos 


(n 1 + to') 


27 rl 


2TV+ 1 


2TV + 1 


(24) 


The Izrailev model in Eqs. (12211241) with a finite number 
of states is considered as the quantum analogue of the 
classical standard mapping on the torus with closed mo¬ 
mentum p and phase 9 , where U nm describes only the 
odd states of the systems, i.e. ip(6) = —tp(—6), pro¬ 
vided we have the case of the quantum resonance, namely 
r = 4irr/(2N + 1), where r is a positive integer, as in 
Eq. ([HI) . The matrix (1241) is obtained by starting the 
derivation from the odd-parity basis of sin(nd) rather 
than the general angular momentum basis exp (in9). 

Nevertheless, we shall use this model for any value of 
r and fc, as a model which in the resonant and in the 
generic case (irrational r/(47r)) corresponds to the clas¬ 
sical kicked rotator, and in the limit TV —> oo approaches 
the infinite dimensional model [Eq. (1 21 [I ] . restricted to 
the symmetry class of the odd eigenfunctions. It is of 


course just one of the possible discrete approximations 
to the continuous infinite dimensional model. 

The difference of behavior between the generic case 
and the quantum resonance shows up only at very large 
times, which grow fast with (2TV+1), as explained above. 
It turns out that also the eigenfunctions and the spectra 
of the eigenphases at finite dimension TV of the matrices 
that we consider do not show any significant differences in 
structural behavior for the rational or irrational t/(47t), 
which we have carefully checked. Indeed, although the 
eigenfunctions and the spectrum of the eigenphases ex¬ 
hibit sensitive dependence on the parameters r and fc, 
their statistical properties are stable against the small 
changes of r and k. This is an advantage, as instead of 
using very large single matrices for the statistical analy¬ 
sis, we can take a large ensemble of smaller matrices for 
values of r and k around some central value of r = To 
and k = ko, which greatly facilitates the numerical cal¬ 
culations and improves the statistical significance of our 
empirical results. Therefore our approach is physically 
mea ning ful. Similar approach was undertaken by Izrailev 
(see (38| and references therein). In Fig. 1 of paper Q 
we show the examples of strongly exponentially local¬ 
ized eigenstates by plotting the natural logarithm of the 
probabilities w n = \u n \ 2 versus the momentum quan¬ 
tum number n, for two different matrix dimensions TV. 
By calculating the localization length £ from the slopes 
cr = 2/C of these eigenfunctions using Eq. (fl5l) we can 
get the first quantitative empirical localization measure 
to be discussed and used later on. The new finding of 
our paper (4ll | is that a has a distribution, which is close 
to the Gaussian (but cannot be exactly that, because a 
is a positive definite quantity). It does not depend on TV 
and survives the limit TV — > oo. Therefore also C has a 
distribution whose variance does not vanish in the limit 
TV —> oo. 

Following dl] and [38] we introduce another measure 
of localization. For each TV-dimensional eigenvector of 
the matrix U nm the information entropy is 

N 

^(« 1 ,..., M W ) = - ^ lnTO„, (25) 

71=1 

where w n = \u n \ 2 , and Y/ n \ u n\ 2 = 1- 

In case of the random matrix theory being applica¬ 
ble to our system [Eqs. (I?T1) and (I22H24I) ]. namely the 
COE (or GOE), due to the isotropic distribution of the 
eigenvectors of a COE of random matrices, we have the 
probability density function of \u n \ on the interval [0,1], 

= 7aW^)75 (1 - K|2,< "" 3>/2 - (26) 

It is easy to show that in the limit TV —> oo this becomes 
a Gaussian distribution 
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and the corresponding information entropy [Eq. (1251) ] is 
equal to 

^n OE = Qw + l) - ^ (0 ^ In QiVa) + Od/N), 

(28) 

where a = exp ( 2 _ 7 ) ~ 0.96, while if; is the digamma 
function and 7 the Euler constant (~ 0.57721...). For 
a uniform distribution over M states w n = 1/M we get 
ft logM, and thus M ft exp (J#n)- Thus, we get the 
insight that the correct measure of localization must be 
proportional to exp(^)v), but properly normalized, such 
that in case of extendedness (GOE/COE) it is equal to 
N. 

Therefore the entropy localization length Ijj is defined 
as 

l H = N exp ( Jf N - Jg$? OE ) . (29) 

Indeed, for entirely extended eigenstates Ih = N. Thus, 
Ih can be calculated for every eigenstate individually. 
However, all eigenstates, while being quite different in 
detail, are exponentially localized, and thus statistically 
very similar. Therefore, in order to minimize the fluctu¬ 
ations one uses the mean localization length d = (Ih), 
which is computed by averaging the entropy over all 
eigenvectors of the same matrix (or even over an ensem¬ 
ble of similar matrices) 


In producing this plot we have used Eqs. (p |) . (p1l for 
A, which is just a rough estimate. Indeed, as we shall 
show below, following 41], A is not a number in a given 
system, at fixed K 1 k and N, but has a distribution, 
whose reciprocal is approximately Gaussian distributed, 
and Eq. msD is just a rough estimate of the mean value 
of A. Therefore we should not be surprised any more to 
see large fluctuations in the scaling law of Fig. [2 an ob¬ 
servation entirely unexplained so far, but clarified in [4l| . 
The statistical properties of the localization measure will 
be discussed below. 

In closing this section we just mention the important 
finding [39] that the level spacing distribution P(S) (of 
the quasienergies) is very well described by the Brody 
distribution fisT 13(1 


Pbk(S) =C 1 S fi exp(-C 2 S fi+1 ) 


(34) 


where the two parameters C\ and C 2 are determined by 
the two generic normalization conditions that must be 
obeyed by any P(S), 

pOO poo 

/ P(S)dS = 1, ( S) = / SP(S)dS = 1, (35) 

Jo Jo 

thus with (S) = 1 being the mean distance between 
neighboring levels (after unfolding). Hence 


Ci = (/3 + 1)C 2 , C 2 = 


r 


P + 2 

P + 1 


1/3+1 


(36) 


d = (Ih) = N exp ((Jif N ) - ^§ OE ) . (30) 

The localization parameter j3\ oc is then defined as 


a _ _ ( 1 h) 

Ploc N ~ N ' 


(31) 


where T(a/) denotes the Gamma function. In the strongly 
localized regime p = 0 we observe Poissonian statistics 
while in the fully chaotic one /? = 1 and the RMT applies. 
The Brody cumulative level spacing distribution is 

Wbr(S) = 1 - exp (~C 2 S p+1 ). (37) 


The parameter that determines the transition from weak 
to strong quantum chaos is neither the strength param¬ 
eter k nor the localization length C , but the ratio of the 
localization length C to the size N of the system in mo¬ 
mentum p 


c _ 1 

N N V t 2 ) 


(32) 


where C ss t\ oc , the theoretical localization length 
Eq. (fl6l) , was derived in ^ • A is the scaling parameter 
of the system. The relationship of A to /?i oc is discussed 
e.g. in [H, m an d further developed in section VII of 
[39f ]. Here we just summarize by showing the empirical 
scaling property of A versus /?i oc in Fig. U where the ap¬ 
proximate analytical description is given by the function 


^i°c = TT—a ’ 7 = 4.04, (33) 

1 + 7 A 

which is similar to the scaling law in in [59:], but not the 
same. Namely, the value 7 = 4.04 differs somewhat from 
7 « 3.2 in (59(, where (3\ oc is plotted versus x ft 4A. 


In Fig. [5] we show the scaling of the Brody spectral pa¬ 
rameter /3 br vs. the localization parameter j3\ oc . They 
are linearly related almost like identity, but we also ob¬ 
serve rather large fluctuations, probably due to the fact 
that /3br is difficult to calculate accurately enough with 
small matrices. 


IV. THE DISTRIBUTION OF THE 
LOCALIZATION MEASURES 

In this section we follow our recent paper [4l| and 
present the results about the distribution of the localiza¬ 
tion measures. Here we restrict our analysis exclusively 
to the case K = 10, as this case is closest to the normal 
diffusion regime fx = 1, as explained in Sec. [TT] First we 
demonstrate that the localization measures C = 2/er and 
Ih are very well defined, linearly related and thus equiv¬ 
alent. In Fig. [T]we show this in the diagram of the mean 
(a) versus 2 /(Ih), where both averagings are over all 
eigenfunctions for matrices of dimension N = 3000, for 7 
nearby values of k around ko, namely k = koztjSk, where 
j = 0,1, 2,3 and 5k = 0.00125, for k 0 = 3,4, 5,..., 19. 
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FIG. 4. [Color online] The parameter /?i oc vs. A for 161 x 398 elements for various values of K and for a wide range of k 
values, where the scaling law [Eq. m] is shown with the black line. 


In the next Fig. [7] we show the relationship of the theo¬ 
retical C in Eq. (1161) and the mean value of the empirical 
2/(a) for k 0 = 3,4, 5,..., 19. It is clearly seen in Fig. [7](a) 
that there are strong fluctuations which we attribute to 
the fact that 2/cr has a certain distribution with nonva¬ 
nishing variance, to be presented and described below, 
and that the theory of C resulting in Eq. © is too sim¬ 
ple, as it corresponds only roughly to the value of 2/ (a). 
On the other hand, in Fig. [T^b) we see again that the 
two empirical localization measures are exactly linearly 
related. We should mention that in the cases of larger 
k > 19 the slopes a are so small, and the localization too 
weak, that we cannot get reliable results. Therefore we 
limit ourselves to the interval 3 < k < 19. 

We have thereby demonstrated that the empirical lo¬ 
calization measures are well defined, while the theoretical 
prediction for their mean values is not good enough. The 
reason is that the localization measures of a given fixed 
system (with fixed K = 10 and k) have a distribution 
with nonvanishing variance, which is out of the scope 
of current semiclassical theories, as they do not predict 
this distribution and the corresponding variance. This 
finding, as the central result of the paper (4l| , is demon¬ 
strated in Fig. [8] The distributions are clearly seen to 
be close to a Gaussian, but cannot be exactly that, as 
a is always a positive definite quantity. Its inverse, the 
localization length equal to C = 2/<r, has a distribution 


whose empirical histograms are much further away from 
a Gaussian, so that in this sense er is the fundamental 
quantity. Indeed, as we will see, it corresponds to the 
finite time Lyapunov exponent known in the theory of 
dynamical systems. 

As Ih and 2/<j are equivalent localization measures, 
the former one is expected also to have a distribution, 
which we demonstrate in the histograms of Fig. [9] 

We have also analyzed how the localization measures 
vary in the semiclassical limit of the increasing value of 
the quantum parameter k , at fixed classical parameter 
K = 10. Indeed, the theoretical estimate of C in Eq. m, 
at fixed K. and remembering k = K/t 1 shows that ap¬ 
proximately the mean value of the localization length 
should increase quadratically with fc, or equivalently, the 
slope a should decrease inversely quadratically with k. 
This prediction is observed, and is demonstrated in the 
Table HI and also in Fig. [TOl 

It is also in agreement with the prediction based on the 
tight-binding approximations in reference (60| [Eq. (6)]. 
We give, in Table Q] the mean slope a and the stan¬ 
dard deviation of er, as well as the mean value of the 
related quantity 2 /Ih and its standard deviation for var¬ 
ious k = ko = 3,4, 5,..., 19, for each of them taking 
seven nearby values of fc, namely k = ko ± jSk , where 
j = 0,1, 2,3 and 6k = 0.00125, for matrices of dimension 
N = 3000. Each histogram for all ko was fitted with the 
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FIG. 5. [Color online] The fit parameter /3br as a function of /?i oc for 161 x 398 elements for various values of K and for a wide 
range of k values. The best fitting straight line is very close to identity. 



FIG. 6. We show (a) versus 2 /(Ih) for matrices of dimension 
N = 3000, for 7 nearby values of k, namely k = ko ± j5k, 
where j = 0,1, 2, 3 and 8k = 0.00125, for ko = 3,4, 5,..., 19. 
The two empirical localization measures are clearly well de¬ 
fined, linearly related and thus equivalent. 


Gaussian distribution and then the mean values and the 
standard deviations were extracted. All four quantities 
decrease to zero with increasing A, meaning that in the 
semiclassical limit the localization lengths monotonically 
increase to infinity, so that in this limit we have asymp¬ 
totically extended states (no localization), and their stan¬ 
dard deviation also goes to zero as 1/A: 2 , which is different 


from the tight-binding approximations in reference |60] 
[Eq. (9)]. 

Next we want to study how does the distribution of 
the localization measure a behave as a function of the 
dimension N of the Izrailev model Eqs. (I22II24I) . Since 
in the limit N —> oo the model converges to the in¬ 
finitely dimensional quantum kicked rotator, we would 
at first sight expect that following the Shepelyansky pic¬ 
ture dil a should converge to its asymptotic value, which 
is sharply defined in the sense that the variance of the 
distribution of a goes to zero inversely with N. Namely, 
at fixed K and A Shepelyansky reduces the problem of 
calculating the localization length to the problem of the 
finite time Lyapunov exponents of the approximate un¬ 
derlying finite dimensional Hamilton system with dimen¬ 
sion 2k. The localization length is then found to be equal 
to the inverse value of the smallest positive Lyapunov ex¬ 
ponent. In our case, the dimension of the matrices N of 
the Izrailev model plays the role of time. As it is known, 
and analyzed in detail in the paper j4l], the finite time 
Lyapunov exponents have a distribution, which is almost 
Gaussian, and its variance decays to zero inversely with 
time. Thus on the basis of this we would expect that the 
variance of a decays inversely with N. 

However, this is not what we observe. In the Table El 
we clearly see that at constant K = 10 and A = 10 the 
mean value of a is constant and obviously equal to its 
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FIG. 7. [Color online] (a) We show £ versus 2 /(a) for matrices 
of dimension N = 1000 (crosses and solid fit line) and for 
matrices of dimension N = 3000 (stars and dashed fit line), for 
7 nearby values of k, namely k = ko ± jSk, where j = 0,1, 2, 3 
and Sk = 0.00125, for ko = 3,4, 5,..., 19. (b) We plot the 
mean value of 2/(A r (cr)) versus /3ioc for ko = 3, 4, 5,..., 19 
and 7 matrices of dimension N = 3000 with k = ko ± jSk, 
where j = 0,1, 2, 3 and the step size Sk = 0.00125. 


asymptotic value of TV = oo, while the variance of a 
does not decrease with N, as 1/JV, but is constant in¬ 
stead, independent of N. This is in disagreement with 
the banded-matrix models of the tight-binding approxi¬ 
mations and thus disagrees with the Eq. (9) of reference 
f60j and also disagrees with the Shepelyansky picture. 
The reason is that the associated Shepelyansky’s Hamil¬ 
ton system is only approximate construction, because 
with increasing N the matrix elements of the Floquet 
propagator (matrix) outside the diagonal band of width 
2k become important, and thus the dimension of the 
Hamilton system cannot be considered finite, constant 
and equal to 2k, but increases with N. As a consequence 
we have the constant value of the variance of <r, and thus 
constant variance of the localization length £ = 2/c, and 
therefore the localization length has a distribution with 
nonvanishing variance even in the limit N = oo. This is 
precisely the reason why the semiclassical prediction of 
the localization length in Eq. m fails in detail and we 
find strong fluctuations in the plot of £ against the 2/cr 
of Figs. S] and 0 The proper theory of the localization 
length must predict its distribution rather than just its 




a 



FIG. 8. We show the histograms of the slopes a for four sys¬ 
tems, matrices of dimension N = 3000, for each of them with 
seven different values of k close to ko = 5, 9,13,17, namely 
k = ko ± jSk, where j = 0, 1,2,3 and Sk = 0.00125: (a) 
ko = 5, (b) ko = 9, (c) ko = 13 and (d) ko = 17. 












13 




FIG. 9. We show the histograms of Ih in (a) and 2/Ih in (b) for the system k = 10 described by the matrices of dimension 
N = 3000. In both cases we show the Gaussian best fit. 



FIG. 10. We show log-log plots in (a) the mean slope (<r) as a function of k, and in (b) the standard deviation of a as a function 
of k. The fitting by a straight line is only on the semiclassical interval 10 < k < 19. In the former case the behavior is roughly 
as 1/k 2 , in agreement with the theoretical estimate 1/k 2 of Eq. (1161 1 . and in the latter case also like 1 /k 2 , surely not as the 
theoretical estimate 1/k based on the Lyapunov exponents method in the reference [60| [Eq. (9)]. 


approximate mean value. 

As it is well known the problem of quantum or dynam¬ 
ical localization is related to the Anderson localization 
model, within the framework of the tight-binding approx¬ 
imation, with hopping transitions between the nearest 
neighbors only. This goes back to the pioneering work 
of Fishman et al. [6lj | . as discussed in Hi , and also 
reviewed in [62}. 

In the paper j4l| we have numerically analyzed the be¬ 
havior of the finite time Lyapunov exponents for a clas¬ 
sical Hamilton system exemplified by the standard map 
(SM), following Fujisaka [63] and Ott [<34| and also for 
the 2x2 random transfer matrices of the tight-binding 
approximation to describe the Anderson localization. In 
both cases we have shown that the distribution of the 
positive Lyapunov exponent is excellently described by a 
Gaussian distribution, whose mean value converges with 
time to the asymptotic value of the infinite time, while 
the variance decays inversely with time t (the number of 
iterations in the case of the standard map), and n, the 
number of random matrices in the product. The latter 
are random unimodular transfer matrices of the tight- 


binding approximation, of the form 

T -( W 1 
" 1 1 0 


(38) 


where W = E — E® is drawn from a distribution, defined 
by a given model. E is the eigenenergy of the system, and 
E® is the fluctuating on-site potential. We have tested 
three quite different distributions for W, namely Gaus¬ 
sian, box distribution and the Cauchy-Lorentz distribu¬ 
tion [41J, and found that the shape of the distribution 


of the positive Lyapunov exponent for any n (=mumber 
of matrices in the product) depends very weakly on the 
overall shape of the W-distribution, while the mean value 
and the variance depend only on the variance of W. In¬ 
deed, the evidence for the predicted decay of the variance 
of the finite time Lyapunov exponents is overwhelming, 
as shown in Fig. 1111 where we plot the standard devia¬ 
tion as a function of time in log-log plot, showing that it 
decays inversely with the square root of time. 

In the context of our Izrailev model the dimension N 
of the matrix plays the role of time. The width of the 
diagonal band is equal to 2k. Shepelyansky reduces the 
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TABLE I. 

The mean value and the standard deviation of the 
slopes a and 2/Ih as a function of 
k = kg = 3, 4, 5,..., 19. For each k = kg we used 


N = 7 x 3000 slopes a (see text). All quantities decay 
to zero in the semiclassical limit,. 


K = 10 - 
k < a > 

N = 7 x 3000 (slopes) - 
SDu < 2 /Ih > 

N = 3000 (2/Z h ) 

sd VIh 

3 

0.06209 

0.01324 

0.062098 

0.01324 

4 

0.04327 

0.01073 

0.043272 

0.01073 

5 

0.04636 

0.00758 

0.046363 

0.00758 

6 

0.04030 

0.00974 

0.040303 

0.00974 

7 

0.04095 

0.00838 

0.040954 

0.00838 

8 

0.03004 

0.00756 

0.030047 

0.00756 

9 

0.03174 

0.00600 

0.031743 

0.00600 

10 

0.02835 

0.00539 

0.028355 

0.00539 

11 

0.02034 

0.00353 

0.020341 

0.00353 

12 

0.02014 

0.00321 

0.020143 

0.00321 

13 

0.01719 

0.0029 

0.017193 

0.00294 

14 

0.01750 

0.00289 

0.017509 

0.00289 

15 

0.01356 

0.00230 

0.013569 

0.00230 

16 

0.01221 

0.00194 

0.012213 

0.00194 

17 

0.00978 

0.00148 

0.009787 

0.00148 

18 

0.00855 

0.00128 

0.008550 

0.00128 

19 

0.00975 

0.00141 

0.009754 

0.00141 


TABLE II. 

The mean value and the variance of the slope a as a 
function of the matrix dimension N for a fixed 
system with K = 10 and k = 10. Both are obviously 
constant. 
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FIG. 11. [Color online] The standard deviation of the positive 
finite time Lyapunov exponents for the standard map (stars) 
and for the product of random transfer matrices with the 
box distribution of W (empty boxes), as a function of time in 
log — log presentation, and their best fits. The slope is exactly 
- 1 / 2 . 


problem of the localization length to the determination 
of the smallest positive Lyapunov exponent (its inverse 
is the localization length) of the underlying finite dimen¬ 
sional Hamilton system with dimension 2k. Then, the 
finite time Lyapunov exponent should have some almost 
Gaussian distribution, whose mean tends to the asymp¬ 
totic Lyapunov exponent with TV —> oo and the variance 
should decrease to zero as 1/N. 

If this picture were exact, then the mean localization 
length as a function of N should converge to the asymp¬ 
totic value, which we do observe in Table lU of Sec. ED 
while the variance does not decay to zero, but rather 
remains constant, independent of N as clearly demon¬ 
strated in Table El From this we conclude that even in 
the limit N —> oo the localization length has a certain dis¬ 
tribution with nonvanishing variance, or more precisely, 
its inverse (the slope a) has an almost Gaussian distribu¬ 
tion with nonvanishing variance. We believe that this is 
the cause of the strong fluctuations observed for example 
in Fig. E and [7fa) of Sec. IIIII 


V. SUMMARY 

We have reviewed our recent results IdoM-flt on the dy¬ 
namical localization in the iV-dimensional Izrailev model. 
The analysis of the classical system (standard map) and 
of the quantum kicked rotator (QKR) has been per¬ 
formed for many different values of the classical kick pa¬ 
rameter K on the interval [K cr . 70] ss [0.97,70], of the 
quantum parameter k and matrix dimensions N in the 
interval [400,3000]. The aspects of classical generally 
anomalous diffusion have been studied and the impor¬ 
tant relevance of the accelerator modes elucidated, and 
the semiclassical approximation for the average localiza¬ 
tion length C has been derived. The entropy localization 
measure Ih has been calculated, and the corresponding 
parameter /3i oc = ( Ih)/N defined. The scaling law of 
Aoc versus A = C/N was found, in agreement with the 
previous results of Izrailev ( [Dsl and references therein). 
However, even after a great numerical effort in maxi¬ 
mally improving the statistical relevance of this scaling 
law, large fluctuations around the mean value have been 
observed. Also, we have shown, that the Brody distri¬ 
bution describes the level spacing distribution very well, 
and the spectral Brody parameter /?br, determining the 
level spacing distribution, was found to be linearly re¬ 
lated to Aoc- 

In the end our main conclusion is the empirical fact 
based on our numerical computations of the eigenfunc¬ 
tions of the iV-dimensional Izrailev model, that the local¬ 
ization length has a distribution with nonvanishing vari¬ 
ance not only for finite N , but even in the limit N —oo. 
This is the reason, we believe, for the strong fluctuations 
in the scaling laws which involve the empirical localiza¬ 
tion measures and the theoretical semiclassical value of 
the localization length. In the Shepelyansky picture [45} 
this might seem to be a contradiction, but the resolution 
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of the puzzle is that in the limit of large N the finite di¬ 
mensional Hamilton system extracted from the Floquet 
propagator of the quantum kicked rotator is not good 
enough, and therefore the matrix elements outside the 
main diagonal band of width 2k play a role, even if they 
are small, but nevertheless plentiful, making the Hamil¬ 
ton system effectively infinite dimensional, with infinitely 
many Lyapunov exponents. This finding is a challenge 
for the improved semiclassical theory of the localization 
length, to derive and explain the discovered distribution 
function. On the other hand, the simple model of the 
Anderson localization based on the tight-binding approx¬ 
imation, with only the nearest neighbor interactions, de¬ 
scribed by the product of 2 x 2 unimodular matrices, has 
a finite dimension, as the transfer matrices are exactly 
two-dimensional, and therefore the variance vanishes in 
the limit of large times n (number of matrices in the 
product) as 1/n. The same conclusion applies to such 
a model with a finite number of interacting neighbors. 
Indeed, according to the references ©[el the variance 
of er should vanish as Var(a) oc 1 /(iVfc 2 ), but our work 
shows that in the quantum kicked rotator this is not ob¬ 


served: the variance does not depend on TV, and decays 
with k faster than 1/fc 2 , namely as 1/k 4 . Thus, we found 
some important differences between the dynamical local¬ 
ization in the quantum kicked rotator and the Anderson 
tight-binding model of localization, and the Shepelyan- 
sky picture, which rest upon the banded matrix models 
with rigorously finite bandwidth. Therefore, the problem 
of calculating the distribution of the localization length 
(or its inverse) in the semiclassical framework is open for 
the future work. Also, the derivation of the Brody dis¬ 
tribution to explain the level spacing distribution of the 
energies [29|, |3l|, 32J in time-independent systems, and of 
the quasienergies [10, III, [39) in time-periodic systems, of 
chaotic eigenstates, is still open for the future. 
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